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ABSTRACT 

An analysis Is made of In-plane and out-of-plane vibra- 
tions of a uniform filament which is bent at right angles 
and the ends of which are fixed. Shear deflection and rota- 
tory inertia are neglected. Curve sheets are presented from 
which one may find the first several frequencies of vibration 
in terms of the geometrical and physical variables. An 
important application is to the vibrations of piping systems. 

Also, as a separate project, a vibration program written 
by G. Pink for the CDC l60^ computer and modified by Y. S. 
Kim, has been rewritten for use with the IBM 360/6? computer. 
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I. INTRODUCTION 



A. GENERAL 

In order to assure long time integrity of nuclear piping, 
the applicable codes, ANSI B31.7 - 1969 Code for Nuclear 
Power Piping and forthcoming edition of Section III of the 
ASME Boiler and Pressure Vessel Code, call for seismic 
analysis of such piping and require that the severity of the 
effects of seismlcally induced vibration be maintained below 
specified acceptable levels. Generally the problem of 
vibration analysis of the major piping systems is so complex 
that a major computational effort is required and this 
demands digital computer analysis using quite sophisticated 
programs. However, in the early stages of a design, when 
alternate layouts are being investigated, it is useful to be 
able to make some reasonable estimates of what a more elab- 
orate analysis would be likely to reveal. 

A recent document [3] concerned with assisting a designer 
to understand and comply with ANSI B31.7 contains material 
summarizing the vibration characteristics of several simple 
piping configurations. It suggests how the designer might 
employ this sort of information to assist him in estimating 
the vibration characteristics of more complicated configur- 
ations. This material should prove to be quite useful. 

However, the scope of the configurations about which 
such information is available is quite limited. One of the 
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basic configurations with which a designer deals is a simple 
ell-shaped configuration. It was decided that it would be 
quite useful to determine and present in a concise form the 
significant results relating to the vibration of such a 
configuration. Specifically, it is assumed that the two 
ends of the configuration are perfectly fixed or clamped, 
that the pipe is of uniform properties along its entire 
length and is of rather small diameter compared to its 
length (hence the word "filament" in the title of this 
thesis) so that shear deflection and rotatory Inertia could 
be neglected, and that the right angle junction of the two 
legs of the ell is small compared to the overall dimensions 
of the ell itself, so that it is not necessary to consider 
the properties of the bend or elbow that would exist in 
practice at the junction. 

Such results are obtained and given herein. The basic 
parameter is the ratio of the length of the shorter leg to 
that of the longer leg; this parameter is designated by the 
Greek letter rho in what follows. A second Important 
parameter, Greek beta in what follows, is required to relate 
the coupling that takes place between bending and torsion in 
the case of what are called out-of-plane vibrations. 

This simple ell-shaped configuration has the great 
virtue of being such that the vibration problem is naturally 
divided into two disjoint sub-parts. One such part involves 
only displacements which lie in the plane of the ell; no 
torsion is involved in these cases. The other part involves 
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both bending and torsion, but the motion of any point on the 
axis of the filament Is only perpendicular to the plane of 
the ell. This gives rise to the two terms used In what 
follows: "In-plane vibrations" and "out-of-plane 

vibrations . " 

B. ASSUMPTIONS AND LIMITATIONS 

It Is well to repeat here specifically the assumptions 
and limitations Involved In this Investigation. 

1. The configuration Is formed from originally straight 
uniform elastic members which meet at a right angle. 

2. The angle at which the two members meet remains a 
right angle throughout the deformation of the 
configuration . 

3. The ends of the configuration are assumed to be 
perfectly clamped. 

The effect of shear deformation has been neglected. 
Only bending and torsional compliances are con- 
sidered and deformations are assumed to be "small." 

5. Rotatory Inertia has been neglected. 

6. Poisson's ratio has been taken to be 0.3. 
Additionally, there are limitations Implicit In the fact 

that values of the roots of the equations which are later 
solved are obtained only to an accuracy of one part In one 
million parts and In the fact that the computational pro- 
grams used a mantissa having 56 bits. These limitations 
resulted In being able to obtain only about the lowest 
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thirteen frequencies for a particular configuration. Also, 
only a few values of beta were Investigated and there are 
some difficulties In the use of linear Interpolation for 
other values of beta. 



C. NOTATION 




A, B, C, D, F, H 


Constants In equation (2-5) and (2-1^1) 


E 


Modulus of elasticity 


G 


Shear modulus of elasticity 




Functions defined In equation (2-6) 


I 


Moment of Inertia of cross section 


J 


Polar moment of Inertia 


L 


Length of pipe 


M 


Bending moment 


P 


Defined In equation (2-^) 


t 


Time 


T 


Torque 


V 


Shear force 


X 


Distance from fixed end 


y 


Lateral deflection 


Y 


Lateral deflection mode shape 


y’ 


Slope 


a 


Defined in equation (2-13) 


B 


Defined in equation (2-29) 


Y 


Phase angle in equation (2-11) 


0 


Angle of twist 




Angle of twist mode shape 
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ij; Constant in equation (2-2) 

y Mass per unit length 

V Mass moment of inertia per unit length 

p Length ratio 

0 ) Circular frequency 
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II. THEORETICAL DEVELOPMENTS 



A. BENDING ANALYSIS 

Neglecting the effects of shear deflection and rotatory 
Inertia the equation governing lateral vibration of a uni- 
form beam Is 




( 2 - 1 ) 



where y is mass per unit length, El is flexural rigidity, 
y is lateral deflection and t is time. Assuming harmonic 
oscillation with circular frequency m, the solution of 
equation (2-1) can be written in the following form: 

y(x,t) = Y(x) cos(ojt + i{;) (2-2) 

where is a constant and Y is a function of x alone. Sub- 
stitution of equation (2-2) into equation (2-1) yields 

_ p^Y(x) = 0 (2-3) 

dx^ 

where 




The general solution of equation (2-3) may be written as 
Y(x) = Ag^ + Bg2 + Cg^ + Dgf^ (2-5) 

where 
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32[Cosh(Px) + Cos(Px)] 



(a) 



St = 



g 2 = 3§[Slnh(Px) - Sln(Px)] 
= ^[Cosh(Px) - Cos(Px)] 
gij = ^[Slnh(Px) + Sin(Px)] 



(b) 

(c) 

(d) 



(2-6) 



At any location x, deflection Y , slope Y , bending moment 
M j and shear force V can be expressed as 

X X 



^x " + °gl^ 

= EIP^CAg^ + Bg|^ + Cg,^ + Dg2) (c) 

= EIP^(Ag^ + Bg^ + Cg2 + Dg^) (d) 

At X = 0 , 



(2-7) 



Yq = A 



; Yq = DP 



(a,b) 



= C(EIP^) ; V. = B(EIp3) (c,d) 



( 2 - 8 ) 



0 " " ’ 0 

Expressing the constants A, B, C and D in terms of Yq, 
Yq, Mq and Vq and substituting into equation (2-7) yields: 



Y 




y’ 




M 




V 


] 

X 



^1 


g^/P 


g3/EIP^ 


g2/EIP 


Pg2 


§1 


gl^/EIP 


g3/EIP‘ 


EIP^g3 


EIPg2 


®1 


gi^/P 


ElP^g^ 


ElP^g^ 


Pg2 


Si 



Y 
M 

V 



(2-9) 
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The above equation expresses the deflection, slope, bending 
moment, and shear force at any location x in terms of those 
at X = 0. 



B. TORSION ANALYSIS 

The equation of motion of a uniform cylinder undergoing 
torsional vibration is. 



where v is mass moment of inertia per unit length, J is 
polar moment of inertia of elastic resistance, and 0 is 
angle of twist. For a sinusoidal solution at circular 
frequency w. 



where y is a constant and cj) Is a function of x alone. Sub- 
stitution In equation (2-10) gives 



0(x,t) = 4>(x) Cos(o)t + y) 



( 2 - 11 ) 




dx 



( 2 - 12 ) 



where 



2 _ VO) 

JG 



2 



a 



(2-13) 



The general solution of equation (2-12) is 



(f)(x) = F Cos(ax) + H Sin(ax) 



( 2 - 1^0 



and torque T Is 



T = = JGal-F Sin(ax) + H Cos (ax)] 



dx 



(2-15) 
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p 



Tq = H(JGa) 



(2-l6a,b) 



4-n = 



Expressing constants P and H in terms of and Tq and 
substituting Into equations (2-l4) and (2-15), the angle of 
twist (f) and torque T at any location x, can be expressed 
in terms of Tq and as follows: 



4) 




Cos (ax) 


Sln(ax)/J Ga 






T 




-JGa Sin(ax) 


Cos (ax) 




T 



(2-17) 



C. IN-PLANE AND OUT-OP-PLANE VIBRATION ANALYSIS OP AN ELL- 
SHAPED UNIPORM PILAMENT WITH CLAMPED ENDS 



1 . In-plane Vibrations 

In Pigure 2-1, an ell-shaped configuration and the 
corresponding free body diagrams are shown. 




Pigure 2-1. Ell-shaped Configuration and Pree 
Body Diagrams for In-plane Vibrations. 



The ratio, p, of the length of the shorter leg to that 



of the longer is a basic geometrical parameter in the 
following analysis. 

At the clamped ends deflections slopes 

YqI, Yq 2 are equal to zero. Also at the junction of the two 
portions, deflections Y^^ and Y ^^2 are equal to zero. Appli- 
cation of equation (2-9) to the two portions yields. 



S: 



EIP 



— M + 
2 01 



EIP 



3 '^01 = 0 



(a) 



EIP 



^^01 + 



g 



EIP 



— V = Y * 

2 '^01 LI 



(b) 



gi| 

g M + —I V = M 
®1 01 p 01 LI 



(c) 



g 



3 .. . 



^ = 0 



2 02 3 02 

EIP ^ EIP^ ^ 



(d) 



gi| gp 

^ ^ = Y, 



Eip2 02 gjp2 02 L2 



(e) 



( 2 - 18 ) 



®1 ^02 p ^02 ^ L2 



(f) 



where g^^^ indicates evaluation of g^^(x) for argument (PL) and 
gj^ indicates evaluation for argument (pPL). 

At the junction. 



LI 



= Y 



L2 



M + M =0 
LI L2 ^ 



Combining equations (2-19) with equations (2-18), 
following set of equations is obtained. 



(2-19a,b) 

the 
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= 0 ( 2 - 20 ) 



S3/eip^ 


0 


g2/EIP^ 


0 




M 

01 


0 


g3/EIP^ 


0 


g2/EIP^ 




M 

‘^02 


S/EIP 


-S4/EIP 


g3/ElP^ 


-g3/EIP^ 




1 — 1 
0 
> 


_Sl 


Si 


64/P 


S4/P 




J02_ 



For a nontrivial solution, the determinant, D, of the above 
matrix must vanish. Elementary manipulations on D lead to 
the following form. 



83 0 S2 0 

0 83 0 g2 

84 -84 S3 -S3 

Si Si 84 84 



0 



(2-21) 



The terms of equation (2-21) are functions of p. For 
any value of p, the zeros (PL)^ of the determinant can be 
found. The corresponding frequencies are evaluated by using 
equation (2-4) as 

“1 (2-22) 

Solution of equation (2-21) will be explained in 
Section III-A. 

2 . Out-of-plane Vibrations 

Out-of-plane vibrations can be analyzed similarly. 
However in this case the effect of torsion must be taken 
into account. 
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In Figure 2-2, free body diagrams of an ell-shaped con- 
figuration for out-of-plane vibrations are shown. 



V. 



LI 




LI 



02 



Figure 2-2. Free Body Diagrams of an Ell-shaped 
Configuration for Out-of-plane Vibrations. 

At the clamped ends deflections slopes 

Yqi, Yq 2 angles of twist 4 >q-|^> 4>q 2 equal to zero. 

Application of equation (2-9) to the two portions yields 
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I 



g- 



+ 






V. 



EIP^ EIP^ 



, _ g, 

'Ll EIP ‘‘01 



Y;. = zrri- + 



V 



EIP 



2 ’01 






*^ L 1 ® 1^01 P ^01 



(a) 



(b) 



(c) 



\l = PgoM^n + gnV 



2 01 ^1 01 



g' 



+ 






V, 



EIP^ EIP^ 



gZi go 

^ V. 



■L2 EIP 02 gjp2 ’02 



^L2 ®1^02 ■*■ P^ ^02 



(d) 

(e ) 

(f) 

(g) 



(2-23) 



\2 = ^§2^02 "■ SiVq2 



(h) 



where indicates evaluation of gj^(x) for argument (PL) and 
g^ indicates evaluation for argument (pPL). 

For the torsion, equation (2-17) applies and gives. 



. _ Sin(aL) ^ 

JGa °1 



= Cos(aL) Tqp 



■ S iU ( up L ) rn 

^L2 JGa ^ 



02 



(a) 

(b) 

(c) 



(2-24) 



Tj^ 2 “ Cos(apL)TQ2 



(d) 



At the junction: 



18 



II 

rH 


• 

C\J 

Fh 


T 

LI 


C\J 

1 

II 


(a,b) 


y' = 

^L1 


*^L2 ’ 


1 — 1 


= y’ 

L2 


(c,d) 


H 

II 


• 

C\J 

> 

\ 


\l 


C\J 

II 


(e,f) 



Specializing to the case of circular cross section and 
assuming Poisson's ratio is 0.3j then 



1 _ 1.3 

JGa Ela 



( 2 - 26 ) 



Combination of equations (2-23), (2-24), (2-25) and 
(2-26) gives the following set of equations. 



0 gijP 

§4? 0 



El 0 



0 

^•^^^ Sin(ocL) 

a 

0 

0 

0 



Cos (aL) 



0 

0 

^^^ -Sln(apL) 

a 

0 

-Cos (apL) 

0 









g2 -g2 




^01 


0 




C\J 

O 


E^ 0 




T 

01 


®1 ®1 




‘''02 






'^Ol 


n 

T ° 




0 ^ 




C\J 

0 

> 


p 







(2-27) 



For a nontrivial solution, the determinant D of this matrix 
must vanish. From equations (2-4), (2-13) and (2-26) it is 
found that 



(aL) = 3(PL)^ 



( 2 - 28 ) 
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where 



6 




1. 3v 






(2-29) 



is a dimensionless parameter. Substituting equation (2-28) 
into equation (2-27) and performing elementary manipulations 
on D give the following equation 



®3 "^3 ° 

0 Sn -^Sin[3(PL)^] 



B 



0 

0 



64 



§2 ®2 



6 - 



i^SintpB(PL)^] g, 0 



= 0 



g 



2 g. 



g]^ 0 



0 

0 



gl gi 



( 2 - 30 ) 



(PL)Cos [Bp(PL)‘^] g^ 0 



0 g^ (PL)Cos [B(PL)^] 



64 



The terms of equation (2-30) are functions of p and B. 
For given values of p and B, the zeros (PL)^ of determinant 
D can be found. The corresponding frequencies are evaluated 
by using equation (2-22). 

Solution of equation (2-30) will be explained in 
Section III-A. 



20 



III. PRESENTATION OF RESULTS 



A. METHODS OF SOLUTION 

The equations (2-21) and (2-30), developed in Chapter II, 
can be solved either by algebraic manipulations or by using 
direct computer solution. 

In a first approach to a solution, algebraic manipula- 
tions were performed to reduce the order of determinants. 

Then transcendental equations were obtained and two computer 
programs, one for the in-plane solution and one for the out- 
of-plane solution, were developed to get the roots of these 
equations. 

In the computer programs the transcendental equations 
were used as function subprograms and roots located in an 
interval. Then the Newton-Raphson method was used to obtain 
good approximate values of the roots. 

The procedure explained above led to some difficulties. 
During the solution of the equation which was derived from 
equation (2-30), the Newton-Raphson method failed to give an 
expected root because the slope of the transcendental function 
reached very high values and this caused a Jump to adjacent 
roots. In one case, the number of iterations required 
exceeded three hundred with an allowable error of 10 (for 
the quantity PL) and increasing the allowable error to 10 ^ 
did not make much difference. Finally, comparison of 
results with the classical solution in the case p is 
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equal to zero showed that some of the roots were 
extraneous . 

The difficulties mentioned in the preceding paragraph 
were eliminated by making direct computer evaluation of 
the determinants (2-21) and (2-30). For this purpose 
subroutine ERODE was developed and each element of the 
determinants was stored in this subprogram. In addition to 
subroutine ERODE, library subroutine DTERM which was 
developed by Jean Eow in 1967 was used to evaluate the 
determinants. This modified procedure resulted in reduction 
of computer time required and in the elimination of extra- 
neous roots. 

Difficulties with occasional failure of the Newton- 
Raphson method were reduced by using the method of successive 
bisectioning. In this method the number of Iterations 
required for one root did not exceed fifty with an allowable 
error of 10”^. This procedure also resulted in conservation 
of computer time. 

Detailed instructions for use of the computer programs 
are given in Appendix E. 

E, GRAPHS 

Numerical results given to a large number of significant 
figures were obtained primarily for the purpose of investi- 
gating the accuracy and Integrity of the solutions. Chapter 
IV will give a discussion of these investigations. 

However, so as to present the results in the most con- 
venient form for practical use, graphical outputs 
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were also obtained by use of the library subroutine 
DRAW. 

One graph was obtained for the in-plane solution. For 
the out-of-plane solution five graphs were obtained cor- 
responding to five different values of B. In each case, 
the ordinates are values of (PL) and abscissas are values 
of p. On these graphs the first (lowest) curve is for the 
first mode, the second curve for the second mode and so 
on. Graphs are given in Figures 3-1 through 3-6. 



C. SAMPLE PROBLEM 



A sample problem will be solved to show how to use the 
graphs given in the preceding article. 

Properties of ell-shaped configuration: 



4 

1-2 

Outside diameter of pipe 
Inside diameter of pipe 
Wall thickness of pipe 
Weight of pipe 
Thickness of insulation 
Weight of insulation 
Weight of contents 
Modulus of elasticity E 

G 



16 * - 6 " 

10 ’ - 3h'' 

6 . 625 " 

6 . 065 " 

0 . 280 " 

18.98 Ib/ft 
1 1/8" 

2.90 Ib/ft 
12.51 Ib/ft 
27,^00,000 Ib/ln^ 
10,600,000 Ib/ln^ 



For the use of graphs, values of p and B must be known. 
For this problem. 
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Figure 3-1. In-plane Vibrations. 
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0.0 0.2 0 . 0.6 0.8 1.0 

Values of Rho 



Figure 3-2. 



Out-of-plane Vibrations, Beta = 0.00001. 
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Values of Rho 



Figure 3-3. Out-of-plane Vibrations, B = 0.001. 
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jure 3-4. 



Out-of-plane Vibrations, 6 



0 . 005 . 
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Figure 3-5. Out-of-plane Vibrations, 6 = 0.01, 
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Values of Rho 



Figure 3-6. Out-of-plane Vibrations, 3 = 0.05. 
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P = = 0.62373737 

^1 



and 3 can be found from equation (2-29). 



3 






T73v 



(3-1) 



where v Is the total mass moment of inertia per unit length 
for the torsion and y is the total mass per unit length 
undergoing lateral vibrations. Assuming that the contents 
(water) do not rotate, then 



Ptotal ^ 1.068875^90 lb-sec^-ln“^ 
and from equation (3-1) 

3 = 0.0150806896 

For evaluation of frequencies, equation (2-22) should be 
used as 



^pipe ^insulation 



^total ^plpe ^insulation ^ '^contents 



numerical values of v 



total ’^total 



Vtotai = 0 . 0509089^105 lb-sec^ 




then 




8.221317961 sec 
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To solve the problem, values of (PL)^^ can be found from 
the graphs. But no graph is available for the value of 3 
found above. Therefore Interpolation is required between 
the values taken from Figure 3-5 and from Figure 3 - 6 . Values 
of (PL)j^ are given in the table below where the exact 
solution was obtained by a computer solution using the 
specified value for 8 . 



MODE 

NO. 


EXACT 

SOLUTION 


SOLUTION OBTAINED BY 
LINEAR INTERPOLATION 


1 


2 . 62^579 


2.624302 


2 


it . 533513 


4.521991 


3 


6.729044 


6.706905 


4 


7 . 862686 


7.817582 


5 


10.518163 


10.230722 


6 


11.908733 


11.698410 


7 


13.730277 


13.442683 


8 


14.259158 


15.619057 


9 


16.376500 


16.379861 


10 


17.185014 


17.403866 



As seen from the table, for the lower modes, linear inter- 
polation gives good results. However there is greater error in using 
linear interpolation for the higher frequencies. One can 
make what appears to be the reasonable conjecture that if 
one were to plot curves with (PL) as ordinate and B as 
abscissa and with different curve sheets corresponding to 
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fixed values of p, there might be some crossing of individual 
curves; that is, for certain values of 6 two distinct modes 
might have the same frequency. Clearly this would confuse 
the significance of characterizing different modes simply 
by ordinal number in the order of Increasing frequency and 
would indicate that a more sophisticated interpolation is 
required. However, no such graphs were prepared and there 
do not seem to be enough values of B represented in the 
graphs at hand to permit making replots. Thus it is 
suggested that some future investigator might find it 
interesting to employ the computer programs, incorporating 
slight modifications, to obtain such graphs and thus inves- 
tigate the reasons for inaccurate interpolation for the 
higher modes. 



I 
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IV. ACCURACY AND ASSURANCE OF SOLUTION 



The allowable error permitted in determining the roots ' 
is 10 ^ and this was assumed to be satisfactory for practical 
purposes. The allowable error can be changed from 10 ”^ to 
any desired value by changing the value EPS in the computer 
programs. This can be done by changing the card 000090 in the 
program IN-PLANE and card 000100 in the program OUT-OF-PLANE. 

Assurance of solution was established in the following ways. 

1 . In the case p is equal to zero, the ell-shaped con- 
figuration becomes a simple built-in beam and solution for 
this case can be compared with the classical solution. 



Degenerate Case; p = 0 


Mode 


Values of (PL) 


No. 


Comparison 
Values * 


In - plane 
Solution 


Out-of-plane 

Solution 


1 


il. 7300^107 


4.7300289 


4.7300411 


2 


7.8532046 


7.8531849 


7.8532047 


3 


10.9956078 


10.9955803 


10.9956081 




14.1371655 


14.1371302 


14.1371655 


5 


17.2787597 


17.2787165 


17.2787594 


6 


20.4203522 


20.4203012 


20.4203526 


7 


23.5619449 


23.5618858 


23.5619449 


8 


26.7035376 


26.7034814 


26.7035412 


9 


29.8451302 


29.8453738 


29.8451861 


10 


32.9867229 


32.9890104 


32.9866973 



* 

Obtained by computer solution of the equation 
Cosh(PL) Cos(PL) = 1 . 
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In equation (2-21) when p Is set to zero the value of 
the determinant Is zero regardless of the value of (PL). 
Therefore p was taken to be 0.00001. 

In out-of-plane case for the larger values of 3j pure 
torsional modes appeared and this can be seen from Figures 
(3-5) and (3-6). The values of (PL) given In the above 
table were taken when 3 Is equal to 0.001. 

2. Consider the cases shown In Figure M-1. 



I 









L 



P = 0.5 

Ca) 



^L 






p = 2.0 



L 






Cb) 



Figure 4-1. Cases Having Reciprocal Values of p. 



The configurations shown In the above figure are identical, 
therefore two systems must have the same frequencies. Using 
equation (2-22) as 



(PL) 



V 



El 

yiT' 



= (PL)2 



i El 



then 



(PL) 

(^ 



al 

bl 



2 
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Mode 


P = 0.5 


0 

CM 

II 

CL 




No. 




(PL)tl 


tplt;7 


1 


^.3755336 


2.1877668 


2.000000000 


2 


7.2490132 


3.6245066 


2.000000000 


3 


8.7635958 


4.3817979 


2.000000000 


4 


10.707676 


5.3538381 


1.999999999 


5 


13.540341 


6.7701705 


2.000000001 


6 


15.024013 


7.5120064 


1.999999999 


7 


16.987974 


8.4939872 


2.000000000 


8 


19.8231918 


9.9115959 


2.000000000 


9 


21.3082077 


10.6541038 


2.000000001 


10 


23.2712829 


11.6356415 


2.000000001 


For out-of 


-plane solution : 


in the case 6 = 


0.001 


Mode 


P = 0.5 


p = 2.0 


(f«al 


No. 




(PL)x^.- 

bx 


(PDbi 


1 


2.9832348 


1.4916187 


1.999998209 


2 


4.7014996 


2.3507534 


1.999996916 


3 


7.3772076 


3.6886158 


1.999993485 


4 


8.8416393 


4.4208958 


1.999965571 


5 


10.8159008 


5 . 4080044 


1.999980037 


6 


13.6035633 


6.8018597 


1.999977062 


7 


15.0968311 


7.5488148 


1.999894233 


8 


17.0625606 


8.5315098 


1.999946210 


9 


19.8673419 


9.9339245 


1.999948965 


10 


21.3568783 


10.6796268 


1.999777575 
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3 . In the case p Is equal to one, it is expected that 



among the in-plane solutions will be many for which the 
moment at the joint vanishes. Thus, these cases also per- 
tain to the clamped-pinned simple beam. Now two clamped- 
pinned simple beams of different lengths, I and Jl*, may have 
matching frequencies in their spectra. That is, for example, 
the second frequency for one may equal the first frequency 
for the other. Specifically the following case was 
considered. 

The third frequency for an equal-leg ell is equal to 
the second frequency of a clamped-pinned beam of length Jl, 
namely 



The first frequency of an equal-leg ell of length SL* is 
equal to the first frequency of clamped-pinned beam of 
length £*, namely 




where 



(PL)^ = 7.06858275 




where 



(PL)^ = 3.92660231 



The ratio P was determined so that then the frequencies are 
equal, thus 




i 



I 



p 



= 0.55550065 



^ (PL), 

a (PL>3 

Then performing an analysis for an ell-shaped configuration 
with the value of p, found above, the second frequency is 
indeed. 




where 

(PD 2 = 7.068384686 

This value shows close agreement with the value 
(PL)^ = 7.0685275 
given above. 
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APPENDIX A 



DESCRIPTION OF PROGRAMS 



A . GENERAL REMARKS 

Programs IN-PLANE and OUT-OF-PLANE are FORTRAN-IV 
language digital computer programs designed for use in the 
vibration analysis of an ell-shaped uniform filament with 
clamped ends . 

Both programs employ double precision arithmetic. Sub- 
routines are provided in binary form. Trigometric and 
hyperbolic functions ^ available in the IBM-360/67 library, 
were used. 



B. PROGRAM STRUCTURE 

Each program consists of a four section main program and 
two subroutines. The function of each section of the main 
program and of each subroutine is given below. 

1. Main program subsections and their functions. 



INPUT Controls read-in of data and allocates 

storage locations for arrays and subscripted 
variables. The allowable error for deter- 
mination of roots is specified in this 
section . 



ITERATION 

AND 

EVALUATION 
OF ROOTS 



In this section the roots of the determinant 
are located in an Interval by searching the 
sign changes of the determinant and then the 
method of successive bisectionlng is used to 
evaluate the roots. The maximum number of 
iterations was limited to three hundred and 
this was assumed sufficient. If the number 
of Iterations becomes as large as three 
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hundred the process terminates and the 
approximate root is tagged to indicate a 
possible Inaccurate result. 


OUTPUT 


Controls output format. 


GRAPH 


Controls the graph output and uses the DRAW 
subroutine provided in the IBM-360/67 
library . 



2. Subroutines and their functions. 



ERODE 


In program IN-PLANE equation (2-21) is stored 
and in program OUT-OF-PLANE equation (2-30) 
is stored. Evaluates each element of the 
determinant and obtains the value of the 
determinant by calling subroutine DTERM. 


DTERM 


This subroutine computes the determinant of 
a matrix of real numbers by Gauss’ method 
of elimination with row pivoting. 


C. PROGRAM 


NOMENCLATURE 


A 


Elements of determinants. 


B 


Checks sign of determinant and if there is a 
sign change moves control to the method of 
successive bisectionlng . 


BB 


Values of 3 in subroutine ERODE. 


BBB 


Checks the acceptability of roots. 


BETA 


Value of 3. 


CEM 


Single precision version of RO. 


DET 


Value of determinant . 


EPS 


Allowable error in roots. 


F 


Value of determinant. 


N 


Index of modes (Ordinal number). 


P 


Values of (PL) used to search Interval. 


R3, R^ 


Values of determinant in the method of 
successive bisectionlng. 


RO 


Values of p. 
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IlS 



4 



ROOT 

RR 

SEV 

X 

XXI ; 
YY 



Roots of determinant (PL). 

Corresponds to RO In subroutine ERODE. 
Single precision version of ROOT. 

Middle of Interval. 

XX2 Upper and lower limits of Interval contain- 
ing the root. 

Values of P, X, XXI and XX2 In subroutine 
ERODE. 
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APPENDIX B 

INSTRUCTIONS FOR USE OF PROGRAMS 
A. DATA REQUIRED 

1. Data required for program IN-PLANE 

a. Number of p*s (NRO). 

b. Number of modes sought (MODE). 

c. Graph control (NGRAPH). 

d. Values of p (RO). 

e. Graph title (ITITLE). 

2 . Data required for program OUT-OF-PLANE 



a. 


Number of 


p's (NRO) 


b. 


Number of 


B's (NBETA). 


c . 


Number of 


modes sought (MODE). 


d. 


Graph control (NGRAPH). 


e . 


Values of 


B (BETA). 


f . 


Values of 


P (RO). 


g- 


Graph title (ITITLE). 



B. DATA DECK ASSEMBLY 
1 . Program IN-PLANE 

The data deck is ordered as follows: 

a. Problem statement card. 

b. Values of p card(s). 

c. Graph title cards. 

Program OUT-OF-PLANE 



ill 



2 . 



The data deck is ordered as follows: 

a. Problem statement card. 

b. Values of 6 card(s). 

c. Values of p card(s). 

d. Graph title cards. 

C. LIMITING THE NUMBER OP ITERATIONS 

In the programs the maximum number of iterations was set 
to 300. If the number of iterations reaches 300 for any 
root, the current value is assumed as a solution but in the 
output it is marked as ITER. EXCEEDED 300. A user either 
accepts the root or changes the cards 000430 in program 
IN-PLANE and 000460 in program OUT-OP-PLANE by changing 300 
to any desired value. 

D. GRAPH 

A user can get graphical output by option. DRAW sub- 
routine, available in IBM- 360/67 library was used and three 
CALL DRAW statements control first, intermediate and last 
curves to be plotted on one graph. DRAW subroutine requires 
single precision arithmetic. Therefore variables CEM and 
SEV are Introduced to change the values of RO and ROOT from 
double precision to single precision. 

E. PORMAT OP DATA 

1. Program IN-PLANE 

a. Problem statement card. 

The card has 3 fields under the control of 



following field specification. 



3 13 



Field names and the Information each field conveys are 
given below In order. 

Field names and content,. 

(1) NRO Number of p*s 

(2) MODE Number of modes sought. 

(3) NGRAPH Code specifying whether there Is 

to be graph output. 

(NGRAPH =0 No graph output) 

(NGRAPH = 1 Graph will be ouiput) 

b. Values of p card(s). 

Field specification 7F10.0. 

If there are more than 7 different values of 
p, additional card(s) will be required with the same field 
specification. 

c. Graph title cards. 

The title of graph is read as data with the 
field specification 6 a 8. For the title two cards are 
required. The first card corresponds to the first line of 
the title and the second card corresponds to the second 
line . 



2 . Program OUT-OF-PLANE 

a. Problem statement card. 

The card has 4 fields under the following field 
specification 

4 13 

Field names are given. below in order. 
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Number of p's. 



(1) NRO 

(2) NBETA Number of 3's. 

(3) MODE Number of modes sought. 

(4) NGRAPH Same as in program IN-PLANE 

b. Values of 3 cards. 

Field specification 7F10.0 and field name is BETA. 

c. Values of p card(s). 

Same as in program IN-PLANE. 

d. Graph title cards. 

Same as in program IN-PLANE. 
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APPENDIX C 



PLOW DIAGRAMS 



A. GENERAL 

Flow diagram of the program OUT-OP-PLANE is given. This 
flow diagram Includes all parts of the main program. Plow 
diagrams of subroutines ERODE and DTERM are not given and 
their functions were explained In Section A-2 In Appendix A. 

Program IN-PLANE Is similar to OUT-OP-PLANE, therefore 
its flow diagram Is not given. 

B. FLOW DIAGRAM OF PROGRAM OUT-OP-PLANE 
Plow diagram follows. 



I 




47 





48 






4 





49 





7 



6 



50 











1 



\ 




-t 



APPENDIX D 



PROGRAM LISTING 



A. INDEX 

1. Program IN-PLANE 
Main program: 



Section 


Sequence 


number 


Page 


INPUT 


000090 




000130 


53 


ITERATION AND 
EVALUATION 
OF ROOTS 


000170 


— 


000550 


53 


OUTPUT 


000590 


- 


000770 


5 ^ 


GRAPH 


000810 


- 


001030 


54 


Subroutines : 










ERODE 


001060 


- 


001300 


55 


DTERM 


001310 


- 


001610 


55 


Program OUT-OF-PLANE 










Main program: 










Section 


Sequence 


number 


Page 


INPUT 


000100 


- 


000150 


57 


ITERATION AND 
EVALUATION 
OF ROOTS 


000190 


— 


000590 


57 


OUTPUT 


000630 




000830 


58 


GRAPH 


000870 




001110 


58 
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I 




ff; 

Pi,- 




Subroutines : 



ERODE 


001140 - 


001600 


59 


DTERM 


001610 - 


001910 


60 



B. LISTING 



© 
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APPENDIX E 



PROGRAM VIPIPE 



A. GENERAL 

As a part of the present thesis investigation, attention 
was given to a digital computer program called VIPIPE which 
was originally written by Fink [1] and later modified and 
augmented by Kim [2] . This program provides solutions for 
in-plane and out-of-plane vibration of a single plane piping 
configuration having quite general geometry (except, of 
course, for its limitation to a single plane). It was 
desired to convert this program from a FORTRAN 63 program 
for use with CDC-l60^ computer formerly employed at the 
Naval Postgraduate School, to a FORTRAN IV program for the 
IBM 360/67 presently in use, in order to provide an indepen- 
dent means of checking the results obtained as indicated in 
the body of this thesis and, also, because it would be 
useful to have a program of such generality available. For 
details of the theory of VIPIPE, refer to the theses written 
by Pink and Kim. 

Because of the limitations on word length in IBM 360/67j 
the converted program, which was called VIPIPE-1, has been 
written to employ double precision arithmetic. Certain 
improvements have been made in the capability to produce 
graphical output. 
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The effort to make this conversion has been found to be 
much greater than was anticipated. It seems that CDC-l 604 
computer is a more "forgiving" machine than is the 
IBM 360/67* In checking out the new program certain bugs 
came to light which were actually present in the earlier 
version but which apparently did not cause trouble. However, 
they represented some syntactical errors that resulted in 
failure on the IBM 36O/67. In some cases it was clear how 
to modify the program so as to eliminate the errors but in 
other cases the difficulty was so deep in the logic that a 
correction could not be made without acquiring as deep a 
knowledge of the program as was possessed by original 
writers. Pink and Kim. 

Thus it cannot be asserted that VIPIPE -1 has been 
successfully converted in every fine detail and it is 
recommended that additional efforts be made to assure the 
overall integrity of the program and to assess its limita- 
tions. It should be noted, however, that no gross errors 
have been found in any results obtained by use of the new 
program; the only question is that of the degree of accuracy 
that it is possible to obtain. 

Comparisons of accuracy are difficult to make. In the 
first place, the word-length is different in the new program, 
and, other things being equal, one would expect the new 
program to provide better accuracy, and to maintain its 
integrity in the face of round-off error beyond the point 
where the earlier programs "bombed out." Thus, agreement 
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between Kim and Fink and disagreement with VIPIPE-1 does not 
mean that the former are "correct" and the latter is in 
error. Furthermore, Kim reports having discovered and cor- 
rected some bugs which he found in Fink’s version. Also, 
there is some difficulty in understanding the precise 
details of some of the problems treated by Fink and Kim. 

With this preface which admonishes the reader not to 
make superficial Inferences from comparisons, in Section C 
some comparisons between results obtained by Kim and results 
obtained using VIPIPE are presented. 

B. DATA DECK ASSEMBLY 

The order of data deck is 

1. Graph title data cards. 

2. Problem statement card. 

3. Boundary condition data card(s). 

Branch joining angle data card(s). 

5. Component extensive properties data card(s). 

6. Component intensive properties data card(s). 

Field specification and format of graph title cards are 

given in Section E-l-c of Appendix B and Information for the 
remaining cards are given in [1] . 

C . SAMPLE PROBLEMS 

Properties of piping systems analyzed in this section 
are given in [1] . Here, only numerical values of frequencies 
will be given. 



1. System 1 (in-plane) 



Mode 

No. 


Comparison* 

Values 


VIPIPE** 
Delta m. (D) 


VIPIPE-1*** 
Delta m. 


1 


75.0999^5^0 


75.09994544 


75.09994535 


2 


207.01589137 


207.01589036 


207.01589099 


3 


^05.83391916 


405.83393183 


405.83391910 


h 


670.86^108383 


670. 86310856 


670.86408365 


5 


1002.15520225 


1002.19978814 


1002.15520252 


6 


1399.70436269 




1399.70434917 


7 


1863.51172599 




1863.51198134 


*Sources of comparison frequencies for system 
are given in Section D. 


1 and system 2 


**Frequencles are given in 


1 [1]. 




*** Frequencies are found 


by VIPIPE-1. 




2. System 2 (In-plane) 






Mode 


Comparison 


VIPIPE 


VIPIPE-1 


No. 


Values 


Delta m. (D) 


Delta m. 


1 


11.80213586 


11.80213389 


11.80213432 


2 


73.96272297 


73.96272167 


73.96272087 


3 


207.09776595 


207.09776685 


207.09776559 




405.82896584 


405.82894132 


405. 82896568 


5 


670.86435913 


670.86200817 


670.86435889 


6 


1002.15518770 


1002.26833752 


1002.15518758 


7 


1399.70436229 


1399.08849534 


1399.70437163 
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3 . System 4 (In-plane) 



Mode 

No. 


VIPIPE 
M-D m. (D) 


VIPIPE -1 
M-D m. 


1 


66.69934902 


67.31642169 


2 


171.22352538 


172.86857155 


3 


335.35007083 


339.09491357 


4 


535.64529414 


544.08929568 


5 


773.44358487 


794.57343086 


4 . System 5 


(In-plane) 




Mode 


VIPIPE 


VIPIPE -1 


No. 


Delta m. (D) 


Delta m. 


1 


116.89479180 


118.95185774 


2 


304.88354398 


302.91288556 


3 


566.43670976 


519.19656352 



Mode 

No. 


VIPIPE 

P.M. m.-C (D) 


VIPIPE -1 
P.M. m. -C 


1 




118.95185779 


2 




302.91288592 


3 




519.19656389 


4 


826.88281929 


794. 49003694 


5 


1126.83797508 


1126.27858174 


6 


1666.58878192 


1464.59530724 


7 


2099.52844030 


1874.97923215 
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D. SOURCES OF COMPARISON FREQUENCIES FOR SYSTEM 1 AND 
SYSTEM 2 



The governing fourth order equation of a homogeneous 
single component straight system 



y 



IV 



El 



y 



(E-1) 



was solved to give 

y = aCosPx + b Cosh Px + c Sin Px + d Sinh Px 

Upon applying boundary conditions for system 1* (y = y* =0 
atx=0, x=L) and simplifying, it was found that the 
eigenvalues of the system must satisfy 

Cos(PL) Cosh(PL) = 1 . (E-2) 

Upon applying the boundary conditions for system 2 (y = y*= 0 
at X = 0, y" = y'" = 0 at x = L), it was found that the 
eigenvalues of the system must satisfy 

Cos(PL) Cosh(PL) = -1 (E-3) 



The roots of equation E-2 and E-3j from which the natural 
frequencies of corresponding systems may easily be found, 
were obtained through Iteration using a digital computer. 



^System 1 is clamped-clamped end beam 
System 2 is clamped-free end beam 



E. PROGRAM LISTING 



Index 

Main program: 
Section 


Sequence 


number 


Pag! 


INPUT • 


000010 




001080 


70 


INVARIANTS 


00II20 




001520 


72 


CONTROL IN PLANE 


001560 




002470 


73 


CONTROL OUT OF PLANE 


002520 




003420 


75 


ITERATION 


003460 




004460 


77 


OUTPUT 


004490 




008135 


79 


Subroutines : 


SUBSEC 


008l40 




008640 


87 


DISTM 


008670 


- 


009220 


88 


DISTMO 


009652 


- 


009676 


89 


RIGID 


009652 


- 


009676 


91 


RIGIO 


009679 


- 


009730 


91 


STIPCO 


009760 




009970 


92 


. STIFOO 


010000 


- 


010230 


92 


STAVEC. 


010260 


- 


010550 


93 


STAVEO 


010580 


- 


010830 


94 


INVERT 


010860 




011570 


94 


HANGER 


OII6OO 


- 


011770 


96 


HANGEO 


011800 


- 


011970 


96 


CPIELO 


012000 


- 


012520 


97 


POINO 


012550 


- 


012770 


98 


CFIELD 


012000 


- 


012520 


98 


POINT 


012550 




013490 


99 
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4 



W 




I 



% 



SPIELD 


013520 


- 


013880 


100 


SPIELO 


013910 


- 


014180 


101 


MATMUL 


014210 




014330 


101 


PINBRA 


014360 




014620 


102 


PINMAT 


014650 


- 


014920 


102 


BRANCH 


014950 


- 


015580 


103 


BRANCO 


015610 




016290 


104 


STATEM 


016320 


- 


016840 


106 


STATEB 


016870 




017610 


107 


BRCOR 


017630 




017940 


108 


DELMA 


017970 


— 


019380 


109 



2 . Listing 
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